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We have acquired a deep i-band image of the BL Lacertae object S5 0716+714 while the target was in an low optical state. Due to the 
faintness of the nucleus, we were able to detect the underlying host galaxy. The host galaxy is measured to have an I-band magnitude 
of 17.5 ± 0.5 and an effective radius of (2.7 ± 0.8) arcsec. Using the host galaxy as a "standard candle", we derive z = 0.31 ± 0.08 
(Icr error) for the host galaxy of S5 0716+714. This redshift is consistent with the redshift z = 0.26 determined by spectroscopy for 3 
galaxies close to S5 0716+714. The effective radius at z = 0.31 would be 12 ± 4 kpc, which is consistent with values obtained for BL 
Lac host galaxies. An optical spectrum acquired during the same epoch shows no identifiable spectral lines. 
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1. Introduction 

BL Lacertae objects (BL Lacs) are active galactic nuclei (AGN) 
characterized by a featureless nonthermal continuum, high opti- 
cal and radio polarization, and variability across the entire elec- 
tromagnetic spectrum. These properties arise from a relativistic 
jet pointed almost towards the observer (Bla ndford & K6nigl| 
|1979). The strong continuum emission from the jet in the opti- 



radio to X-ray bands ( Wagner et al.|1996l|Ostorero et al.|2006| l. 
The nucleus of S5 0716+714 is typically bright in the optical 



cal sometimes presents a problem when attempting to properly 
characterize the underlying host galaxy or measure its redshift, 
especially since by definition the emission lines are very weak 
in BL Lacs. Hence, it is unsurprising that even after very deep 
searches (e.g. |Sbarufatti et al.|2006j ) many BL Lacs remain with- 
out a spectroscopically determined redshift. 

The redshift is a fundamental property of any extragalactic 
object, without which e.g. the luminosity or intrinsic variabil- 
ity properties of the target cannot be determined. Furthermore, 
many BL Lacs have been recently detected at TeV gamma-rays 
using ground-based Cherenkov telescopes MAGIC, HESS, and 
VERITAS (e.g. Hinton 2007 1. Since VHE gamma-rays can be 
absorbed by the interaction with low energy photons of the EBL 
via pair production, this opens the possibility of measuring the 
amount of extragalactic background light (EBL, |Nikishov|1962[ 
|Stecker et al.||1992^ ) in the optical through to the far infrared, 
which provides important information about the galaxy and star 
formation history. The absorption depends strongly on the dis- 
tance of the source and the energy of the gamma-rays. If the 
redshift of the source is known, the VHE gamma-ray spectrum 
can be used to derive limits on the EBL (e.g. jAharonian et al. 
[20061 IMazin & Rauel[2007l l. On the other hand, if the redshift 
is unknown, the VHE gamma-ray spectra can be used to set an 
upper limit to the redshift of the source (e.g. Mazin & Goebel 

After its discovery as a bright (Sschz > 1 Jy) flat-spectrum 
(a < -0.5, Sv oc v") radio source, the BL Lac object S5 
0716+714 has been studied intensively at all frequency bands. 
The object is highly variable with rapid variations observed from 



(see e.g.jNescietal 



acterize its host ga 



Pursimo et al. 2002 
icaUy dStickel et al. 



2005 1, thus previous attempts to either char- 
axy ( fstickel e"raL][T993) [Urry et al.|[2000l 



I or to determine its redshift spectroscop- 



1993[ [Marcha et aL||1996^ have not been 



successful. Thus, in spite of numerous variability studies of S5 
0716+714, it has never been possible to determine reliably e.g. 
the linear dimensions and luminosities of the varying compo- 
nents. 

S5 0716+714 has also been detected in TeV gamma-rays 
(|Teshima et al.| |2008). Given that previous estimates (e.g. 



'Wa gner et al.||1996 ) place S5 0716+714 at z > 0.3, this would 
make the target one of the most distant TeV sources detected so 
far. Given the obvious importance of S5 0716+714 to EBL stud- 
ies, it is be important to secure an accurate spectroscopic redshift 
for S5 0716+714 or at least constrain the redshift by some other 
means. 

S5 0716+714 is regularly monitored at Tuorla Observatory 
as part of the blazar monitoring progranr] On December 17, 
2007 we observed S5 0716+714 to go into a fairly deep opti- 
cal minimum (R ~ 14.8) and initiated prompt i-band imaging 
at the Nordic Optical Telescope (NOT) to detect its host galaxy. 
The imaging was performed five days after the minimum and the 
results are reported in this Letter. 

Throughout this Letter, we use the cosmology Ho - 70 km 
s"' Mpc"', Qm = 0.3 and Qa = 0.7. 



2. Observations and data analysis 

S5 0716+714 was observed at the Nordic Optical Telescope 
(NOT), La Palma on Dec 22, 2007. We used the ALFOSC in- 
strument in imaging mode with a 2k E2V CCD chip of a gain 
0.726 e" ADU"' and readout noise 3.2 electrons. Twenty-eight 
exposures of 45 s were acquired using an i-band interference fil- 



http://users.utu.fi/~kani/lm/index.html 
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Fig. 1. The summed i-band image of S5 0716+714 obtained at 
the NOT in Dec 22, 2007. The field size is 3.85 x 3.7 arcmin, 
north is up and east is to the left. Object T is a transient or a very 
red object not visible in previous images of S5 0716+714, most 
of which were obtained in the R-band. The other labeled objects 
are discussed in the text. 



ter of almost uniform transmission between 725 and 875 nm, 
which provided a total exposure time of 1260 s. The i-band filter 
was chosen because it detects light above the 4000 A break of 
the host galaxy even at z = 0.5, which maximizes the amount of 
light from the host galaxy. 

The individual images were bias-subtracted and flat-fielded 
with twilight flats in the usual way using IRAFrl A fringe pat- 
tern with an amplitude of ~ 2.5% of the sky background was 
visible in the reduced images. A fringe correction image was 
produced by median combining all 28 images and applying si- 
multaneously a sigma clipping cut to the pixel intensity values. 
Since the position of S5 0716+714 on the CCD was altered be- 
tween the exposures, this procedure produced a clean correction 
image consisting purely of the fringe pattern. After correction 
for the fringe pattern, the individual frames were registered and 
summed. The summed image shown in Figure [T] has a FWHM 
of 0.80 arcsec and the sky background is flat within 0.4% over 
the entire field of view. 

The field was calibrated using stars 3, 5, and 6 for which 
Ghisellini et al. ( |1997| ) published Cousins 1-band magnitudes. 
Unfortunately, all three stars were saturated in their cores (r 
< (jf.'A - 0'.'8), and could not be used directly for calibration 
purposes. Instead, we used a 45 s i-band exposure of the field 
acquired immediately before the S5 0716+714 sequence. The 
FWHM in this image is significantly worse than in the the S5 
0716+714 sequence and stars 3, 5, and 6 are not saturated. Using 
these stars we first determined the zero point of the 45 s expo- 
sure. All three stars gave consistent results, the maximum dif- 
ference between the zero points was 0.047 mag, consistent with 
the errors of the Ic magnitudes 0.04-0.05 mag in Ghisellini et al. 



^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



( |1997| l. Using this zero point, we determined the Ic magnitudes 
of stars SI, S2, and S3 and from these three stars the zero point 
of the combined long exposure. The uncertainly in the zero point 
was dominated by the uncertainty in the Ic magnitudes of stars 3, 
5, and 6 and by the color effects between our i-band filter and the 
Ic filter We estimate the former contribution to the uncertainty 
to be 0.03 mag. 

To estimate the magnitude of the color effects, we used the 
established relation between Ic and and SDSS i given by Lupton 
(2005) in the SDSS web pages [] / = / - 0.3780(/ - z) - 0.3974. 
The SDSS i has approximately the same width as our i-band 
filter, but the central wavelength is ~ 50 nm lower For stars 3, 5, 
and 6, the i-z color is ~ -0.01 - 0.06. Assuming stars SI -S3 have 
similar colors and since an elliptical galaxy at z = 0.3-0.5 has a 
color i - z ~ 0.4 - 0.5, the color effects between the calibration 
stars and the host galaxy of S5 0716+714 may cause an error 
of up to 0.2 mag in SDSS i. Since our filter is close to SDSS 
i, we expect the color effects to be similar in value. Since the 
color errors dominate the total eiTor of the zero point, we assign 
a formal error of 0.2 mag to the magnitude zero point. 

The host galaxy was analyzed by fitting a two-dimensional 
surface brightness model to the light distribution of S5 
0716+714. Details of this process can be found in INilsson 
|et al. (1999 1. In short, the model consists of two components, 
an unresolved core representing the BL Lac nucleus and a de 
Vaucouleurs profile representing the host galaxy. Previous imag- 
ing (e.g. Urry et al. 2000 1 has shown that the de Vaucouleurs pro- 



file describes well the surface brightness profiles of BL Lac host 
galaxies, although deviations from this law have also been ob- 
served (e.g. Nilsson et al. 2003| l. For our purposes, it is sufficient 
to know that BL Lac host galaxies are bulge-dominated systems 
and no BL Lacertae object has ever been reliably associated with 
a disk type host. 

The model is described by 5 parameters: position (x,y), core 
magnitude nic, host galaxy magnitude nin, and host galaxy ef- 
fective radius reff. The model parameters are adjusted using an 
iterative Levenberg-Marquardt loop, until the minimum value of 
the x^ statistic between the model and the data is found. Only 
pixels within 9'.'5 from the center of S5 0716+714 were included 
in the fit and pixels influenced by an I = 20.8 mag nearby com- 
panion 4'.'2 west of S5 0716+714 were excluded. 

The model is convolved by the observed PSF, determined 
from stars in the vicinity of S5 0716+714. Since stars 5 and 6 
were saturated in their cores (r < 0.4 arcsec; the peak counts of 
S5 0716+714 remained under 58 000 ADUs in afl images.), we 
constructed the PSF from stars 5 and SI by joining smoothly 
the outer parts of star 5 with the inner part of S 1 (hereafter we 
refer to this PSF as PSFl). A second PSF was constructed in a 
similar way from stars 6 and S2 (PSF2) to study the effect of 
PSF variability over the field of view. 

Since the host galaxy is faint compared to the nucleus, the 
results are highly dependent on the accuracy of the PSF and its 
stability over the field of view. We analyzed the PSF stability by 
extracting the surface brightness profiles of stars 2, 3, 5, 6, SI, 
and S2 and scaling them to the same magnitude. In the lower 
panel of Fig. [2] we show the surface brightness profiles of these 
stars relative to star 5, which is equal to PSFl at r > 0'.'4. We note 
that in this representation star 5 appears as a horizontal line. The 
surface brightness profiles of the stars agree to within 0.2 mag 
all the way to the outer fitting radius of 9'.'5. The rms scatter be- 
tween the surface brightness profiles is 0.8% close to the center 



^ http://www.sdss.org/dr4/algorithms/sdssUBVRITransform. 
html#Lupton2005 
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Table 1. Model fit results. PSFl was derived from stars 5 and 
SI in Fig.[T]and PSF2 from stars 6 and S2. 
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and increases to ~ 10% at the outer fitting radius. Part of this 
increase is undoubtedly due to the increase in noise in measur- 
ing the surface brightness at the faint outer wings of the stars. 
However, we assume conservatively that the increase is entirely 
due to PSF variability and represent it by a parabolic expression 

o-psF(r) = 0.008 + 4 X 10"^^ ■ r^ , 

where crps,p{r) is the uncertainty in the PSF relative to the inten- 
sity at radius r (pixels) from the PSF center. We then compute 
the expected variance of a pixel cr- with intensity / (ADU) and 
distance r from the center of S5 0716-1-714 from the expression 



(^2 = 






+ [o-psF(r) * if 



where G is the effective gain and R is the effective readout noise, 
and compute the x^ from 



X' 



-Z 



(li - Mif 



where M is the model intensity and the summation is over all 
unmasked pixels within the fitting radius. 

In addition to the i-band images, two spectra with an expo- 
sure time of 900 s were obtained using grism #4 of ALFOSC, 
which covers the wavelength range 3200-9100 A. The spectral 
resolution achieved using the r.'3 slit was 21 A. Due to the effect 
of the secondary spectrum and fringing longwards of 6000 A, we 
studied only the wavelength range 4000-6000 A. The summed 
spectrum has a signal-to-noise ratio (S/N) of ~ 250, but the spec- 
trum remains featureless. 

3. Results 

Tablefllindicates the results of the model fits, and the upper panel 
of Fig^shows the one-dimensional surface brightness profile of 
S5 07164-714 and the model profiles. The I-band magnitude of 



the nucleus translates into R ~ 14.2 using (R-I) = 0.5 ( |Gu et al. 
[2006,, Stalin et al.|2006i , so the nucleus had brightened by ~ 0.6 
mag in the R-band from the deep minimum on Dec 17th. 

The host galaxy is faint (I - 17.5) compared to the nucleus 
and never exceeds the core brightness at any radius, but is vis- 
ible as a small excess in the surface brightness profile of S5 
0716H-714 (see Fig.l2|. Since PSFl was constructed from stars 
closer to S5 0716H-714 than PSF2, it is unsurprising that PSFl 
provides a more accurate description of the data (smaller y^) 
than PSF2. The fit with core + galaxy reduces the x^ of the fit, 
but none of the fits provide formally satisfactory fits, most likely 
due to higher-order PSF variability unaccounted for in our radial 
PSF variability model. 

We completed several tests to ensure that the observed ex- 
cess is real and not due to PSF variability over the field of view. 
Firstly, we compared the scaled surface brightness profiles of 
stars 2, 3, 6, SI, and S2 and S5 0716-1-714 with the surface 
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Fig. 2. Upper panel: The I-band surface brightness profile of S5 
0716-1-714 (filled symbols). The solid line shows the core + host 
galaxy model, the dashed line the core component and the dot- 
dashed line the host galaxy component. The vertical dashed line 
shows the outer radius of the model-fitting region. Lower panel: 
The surface brightness profiles of S5 0716+714 (symbols and 
solid line) and stars 2, 3, 5, 6, SI, and S2 (dashed lines) relative 
to star 5. Each star is traced out to a radius at which the surface 
brightness can be determined to higher accuracy than 10%. 



brightness profile of star 5 (see the lower panel of Fig. 2). There 
is a clear excess around S5 0716+714 that exceeds the rms scat- 
ter of stellar profiles by a factor of 3-6 at r > 2". The excess is 
clearly above any PSF variability observed over the field of view. 

Secondly, although the fits with PSF2 provide a slightly 
worse fit than with PSFl, the host galaxy component is still 
present with similar brightness and effective radius as with PSFl . 
Also this test indicates that PSF variability does not affect the 
results significantly. Finally, we fitted core + host galaxy mod- 
els to stars 6, S2, and S4, which are similar in brightness to 
S5 0716+714. Using PSFl, the fit converged towards a solution 
'"eff — > in all three cases, i.e. no host galaxy component was 
detected in targets known to be unresolved. Based on the three 
test mentioned above, we conclude that the excess around S5 
0716+714 is real and we have detected the host galaxy. 

Knowing the host galaxy magnitude, we can estimate the 
redshift of S5 0716+714, albeit with rather large margins. 
ISbarufatti et al.|p005] l demonstrated that the distribution of BL 
Lac host galaxy absolute magnitude Mr is almost Gaussian with 
an average of -22.8 and cr - 0.5, using the same cosmology 
as we do here, and BL Lac host galaxies can therefore be used 
as a "standard candle" to estimate distances. Before using their 
method, we first have to transform the apparent I-band magni- 
tude of the host galaxy to apparent R-band magnitude. We cor- 
rect the host galaxy magnitude I = 17.47 for the galactic absorp- 
tion by assuming A/ = 0.06 (Schlegel et al. [T998| l. Since the 
observed R-I color depends on redshift, we have to determine 
the redshift by iteration, starting from z = and using Eq. (2) in 
ISbarufatti et al.|p005| l and typical elliptical galaxy colors from 
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[Fukugita et"aLl ( |1995| ). This iteration yields R = 18.3 and z = 
0.31 in the adopted cosmology. 

To compute the uncertainty in the derived redshift, we use 
the difference between the fits with PSFl and PSF2 as an esti- 
mate for the Icr fitting uncertainty of the host galaxy magnitude, 
i.e. a-j i,„si = 0.4 mag. Summing this in quadrature with the error 
of the zero point (0.2 mag), the total uncertainty in the I-band 
magnitude of the host is 0.45 mag. Performing the redshift iter- 
ation at /host + 0.45 mag and /host - 0.45 mag yields an error of 
+0.5 mag in the R-band magnitude of the host and ±0.06 for z. 
By adding to this in quadrature the inherent uncertainty in the 
method (Az = 0.05, [Sbarufatti et al.||2005[ ), we derive the final 
error in z to be 0.08. 

Finally, we mention that the effective radius of the host 
galaxy 2'.'7 ± 0'.'9 translates into 12+4 kpc at z = 0.31, which 
is consistent with typical values found for BL Lac host galaxies. 

4. Discussion 

Our measurement of redshift does not significantly differ from 
the value z - 0.26 obtained for three galaxies close to S5 
0716+714 using spectroscopy (Stickel et al. 1993[ Bychkova 
et al. 2006 galaxies G1-G3 in Fig. IIJ. Both the host galaxy 
and nearby environment therefore have properties indicative of 
an object at z = 0.26. The nearby companion 4'.'2 W of S5 
0716+714 remains unresolved in our image. 

There have been four previous attempts to derive the host 
galaxy properties of S5 0716+714 with varying success. IStickel 
et al.| ( |1993| l found the object to be unresolved in a 500 s R- 
band image obtained under 1.1 arcsec seeing at the Calar Alto 



3.5 m telescope. , Scarpa et al. (2000l obtained a 614 s HST im- 
age through the F702W filter and found the surface brightness 
profile to be consistent with the PSF and derived an upper limit 
R > 20.0 for the host galaxy. Pursimo et al. ( [2002 ) found that 
S5 0716+714 was unresolved in a 2120 s R-band exposure at 
the NOT, but the PSF was determined partly from a different im- 
age and the PSF shape was uncertain. Finally, [Bychkova et al.| 
( [200 6') reported a detection of excess radiation in the wings of 
S5 0716+714 in VRI images of exposure times between 600 and 
1500 s obtained with the 6 m SAO telescope. They described the 
excess as very large (15"x 7") and elliptical in shape. 

Bychkova et al. ( 2006[l, hereafter B06, did not perform model 



fitting to their image but rather derived the host galaxy proper- 
ties from PSF-subtracted images. Since neither surface bright- 
ness values nor the FWHM were provided by B06, it is difficult 
to assess if the large extent of the host galaxy reported in B06 is 
consistent with our host galaxy magnitude and effective radius. 
However, they indicated that I ~ 15.1 for the extended emission, 
far too bright to be consistent with the faint host galaxy that we 
observe in our i-band image. 



The upper limit R > 20.0 derived by Scarpa et al. (2000i, 
hereafter SOO, may at first glance appear to contradict our es- 
timate for the host galaxy of R = (18.3 + 0.5). Based on this 
upper limit, ,Sbai-ufatti et al.| ( |2005| l infeiTed that z > 0.52 for S5 
0716+714. The brightness of the nucleus of S5 0716+714 was 
the same in SOO (R = 14.18) as here (R ~ 14.2), so the differ- 
ence cannot be explained by a greater core dominance in one of 
our studies. However, the upper limits in |Scarpa et al.| ( [2000| ) are 
based on statistical errors and do not take into account possible 
systematic errors. Consequently, there might be considerable un- 
certainty in the derived upper limit and the two results could be 
consistent within errors. 

We finally note that our redshift estimate is based on the as- 
sumption that the host galaxy of S5 0716+714 has an average lu- 



minosity. The failed previous attempts to detect the host galaxy 
suggest that the host galaxy may actually be below average lu- 
minosity and the redshift consequently lower that 0.31. We note 
that the scatter in the BL Lac host galaxy luminosities is already 
included in our redshift error since this scatter is included in the 
0.05 uncertainty of the method by Sbarufatti et al. ( 2005[ l. Thus, 
within a margin of error of 2cr, the redshift could be in the range 
of0.15to0.47. 



5. Conclusions 

A deep i-band image of the BL Lacertae object S5 0716+714 
acquired at the Nordic Optical Telescope while the target was 
in an low optical state has enabled us to detect the underlying 
host galaxy. The host galaxy has an I-band magnitude of 17.5 
+ 0.5 and an effective radius of (2.7 ± 0.8) arcsec. Using the 
host galaxy as a "standard candle" and a method proposed by 
[Sbarufatti et "all 120051 ), we derive a redshift of z = 0.31 ± 0.08 
(1(T error) for the host galaxy of S5 0716+714. This redshift is 
consistent with the redshift z = 0.26 determined by spectroscopy 
for 3 galaxies close to S5 0716+714. The corresponding effec- 
tive radius of the host galaxy at z = 0.31 would be 12 + 4 kpc, 
also consistent with the values obtained for BL Lac host galaxies 
in other studies. An optical spectrum obtained at the same epoch 
shows a featureless continuum with no identifiable spectral lines. 
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